Chemical engineering of adamantane by lithium functionalization: A first-principles 

density functional theory study 
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Using first-principle density functional theory, we investigated the hydrogen storage capacity of Li 
functionalized adamantane. We showed that if one of the acidic hydrogen atoms of adamantane is 
replaced by Li/Li + , the resulting complex is activated and ready to adsorb hydrogen molecules at a 
high gravimetric weight percent of around ~7.0 %. Due to polarization of hydrogen molecules under 
the induced electric field generated by positively charged Li/Li + , they are adsorbed on ADM.Li/Li + 
complexes with an average binding energy of ~-0.15 eV/fi2, desirable for hydrogen storage appli- 
cations. We also examined the possibility of the replacement of a larger number of acidic hydrogen 
atoms of adamantane by Li/Li + and the possibility of aggregations of formed complexes in ex- 
periments. The stabilities of the proposed structures were investigated by calculating vibrational 
spectra and doing MD simulations. 

PACS numbers: 88.30. R-,81.05.U-,81.07.Nb,87.85.Qr 



I. INTRODUCTION 

Hydrogen storage is of great interest as environmen- 
tally clean and efficient fuels required for future energy 
applications. Several pioneering strategies have been de- 
veloped and significant performances have been achieved 
for hydrogen storage, including chemisorption of dihy- 
drogen in the form of light metal hydrides, metal nitrides 
and imides, and physisorption of dihydrogen onto carbon, 
clathrate hydrates, and porous network materials such as 
carbon nanotubes, zeolites, and metal-organic framework 
(MOF) materials 1 . However, hydrogen storage in these 
systems requires either high pressure or very low temper- 
ature, or both, thus severely limiting the applicability for 
mobile applications. Furthermore, their synthesis at bulk 
scale in a discrete and tailored fashion in a high yield is 
quite difficult 2 ^. Thus, the synthesis of functional mate- 
rials with high hydrogen uptake and delivery under safe 
and ambient conditions remain a key challenge for estab- 
lishing a hydrogen economy. 

To improve the capability of hydrogen storage in ma- 
terials, it has been suggested that they should be doped 
with transition metal impurity or alkali metals. Gener- 
ally, however, metal impurities bind to carbon surfaces 
weakly and they undergo aggregation after several sub- 
sequent fueling cycles^. Also synthesis of lithium doped 
materials leads to the agglomeration of Li atoms, result- 
ing in uneven binding of Li on the metal surface, thereby 
reducing the storage capacity of the materials. To date, 
no material that consists of high Li content with ul- 
tra high hydrogen storage capacity has been reported^. 
Hence, theoretical suggestions and speculations thus far 
proposed have not yielded experimental or technologi- 
cal methods for large-scale production. Diamondoids are 
hydrocarbons with cubic-diamond cage structures that 



offer valuable chemical building blocks for new materi- 
als^. Since adamantane (ADM) can be isolated in large 
amounts from crude oil, research activities focusing on 
such compounds as novel materials are rapidly emerging. 
Experimental and theoretical studies on ADM showed 
that these molecules undergo self-assembly on metal sur- 
faces or into molecular crystals in vacuum with high 
porosity 9 - - — . Furthermore, a synthetic approach to func- 
tionalize the ADM molecule with Li in high purity and 
yield can be realized by the addition of bases such as 
organolithiumi^ and the product has been well charac- 
terized using various spectroscopic methods 18 . 

In this paper we have considered the hydrogen stor- 
age application of Li functionalized ADM. Our first- 
principles calculations show that by chemical substitu- 
tions of hydrogen atoms of ADM with Li/Li + , it becomes 
activated and ready to adsorb hydrogen molecules with 
binding energy of -0.1 to -0.2 eV, with high gravimetric 
weight percents of more than 7.0 %. The Department Of 
Energy (DOE) targets for economical hydrogen storage 
materials specified for 2010 are 6.0 wt % and 45 kg/m 3 for 
gravimetric and volumetric capacities, respectivel y 19 ! 20 . 
As another important criterion, recently, Bhatia and My- 
ers studied the optimum thermodynamic conditions for 
hydrogen adsorption, employing the Langmuir equation 
and derived relationships between the operating pressure 
of a storage tank and the enthalpy of adsorption required 
for storage near room temperature. They have found that 
the average optimal adsorption enthalpy should be in the 



range of 0.1-0.2 eV/H 



21,22 



It is seen that our designed 



materials satisfy the criteria for the gravimetric density of 
hydrogen in storage media by the DOE as well as above 
predicted optimal energy window for hydrogen adsorp- 
tion. Therefore diamondoids-Li/Li + complexes might be 
good candidates for hydrogen storage materials. 
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II. COMPUTATIONAL METHOD 

Density functional theory and ab-initio molecular dy- 
namics calculations were performed with the Vienna ab- 
initio Simulation Package (VASP)^ 3 - using the projec- 
tor augmented wave (PAW) method^ to describe the 
ion-electron interactions. Electron exchange-correlation 
functionals were represented with the generalized gra- 
dient approximation (GGA), with two nonlocal correc- 
tions, Perdew and Wang (PW91)2£ and Perdew-Burke- 
Ernzerhof (PBE)? ? functionals. We used a simple large 
cubic cell with a lattice parameter of 35 A. For all calcula- 
tions, a plane-wave expansion cut-off of 400 eV was used 
and the surface Brillouin zone integration was calculated 
using a T point due to the use of large cubic cell size and 
cluster-type calculation. The structures were fully opti- 
mized until the magnitude of force on each ion becames 
less than 0.005 eV/A . The convergence criterion on the 
total energy was set to lxlO -5 eV. Molecular dynamics 
calculations were carried out at a constant temperature 
for 48 ps with a time scale of 2 fs to determine the sta- 
bility of the alkali functionalized ADM molecules. The 
calculations with hybrid meta GGA functional, (M05- 
2X)22, and the second-order M011er-Plesset perturbation 
theory (MP2}2& were carried out as implemented in the 
GAUSSIAN09 program package^. DFT and MP 2 cal- 
culations were carried out using the 6-31+G(d, p) and 
6-311+G(2df, p) basis sets^o. The choice of M05-2X/6- 
311+G(2df, p) for DFT calculations is justified as a com- 
promise between reliable results and a reasonable com- 
putational cost compared to MP2 method. 

The consecutive binding energy of Li, the required en- 
ergy to separate Li from ADM.Li m complex, is calculated 
subsequently as E b = £ADM.Li m+1 - £^ADM.Li m - E u , 
where -Eadm.ll„ + i, ^ADM.Li m {m > 1) are the total en- 
ergies of ADM with m+1 and m substituted Li atoms, 
respectively. Eu is the total energy of a Li atom. The 
binding energy of Li in the ADM. Li is calculated as 
Eb = -Eadm.li~-Ea.dm* —Eu, where -Eadm.li and -Eadm* 
are the total energy of optimized structures of ADM. Li 
and ADM with one detached acidic hydrogen atom, re- 
spectively. Similarly, the binding energy of the Li + in 
ADM.Li + is calculated as Eb = E ABM Li + — -Eadm* — 

The binding energy per hydrogen molecule, Eb, was 
calculated using 

Eb = - [-EADM.Li m -(H 2 )„ _ -EADM.Li m - «-Eh 2 ] , (1) 

where £ADM.Li ro -(H 2 ) n , E A DM.u m and En 2 are the total 
energies obtained for alkali atom-doped ADM containing 
n hydrogen molecule(s) and m lithium atom(s), the alkali 
atom-doped ADM system and an isolated H2 molecule 
being located in the same supercell. 

In our study, the excess and depletion charge or 
difference charge density Ap is estimated by Ap = 
p(ADM.Li)-p(ADM*)-p(Li), where p(ADM.Li) stands 
for the charge density of full relaxed structure of ADM. Li. 



p(ADM*) and p(Li) are the charge densities of ADM* 
(ADM. Li with detached Li) and Li atom, which are ob- 
tained from two separate single point energy calcula- 
tions without any relaxation while the positions of all the 
atoms in ADM* and Li are kept fixed as the positions of 
their corresponding atoms in ADM. Li. 



III. RESULTS AND DISCUSSION 

Before evaluating the hydrogen storage property of 
chemically modified diamondoid-based complexes, as an 
example of pristine diamondoid structures, we exam- 
ined the hydrogen storage property of ADM using var- 
ious density functional approaches. In all our calcula- 
tions, it was observed that due to the very high stabil- 
ity of ADM 31 , the binding energy of hydrogen molecules 
on such compounds was very small, i.e., on the or- 
der of ~ —0.001 eV, which is out of the energy range, 
from -0.1 to -0.2 eV, desirable for reversible H2 ad- 
sorption/desorption near room temperature for hydro- 
gen storage applications^^. Hence, pristine diamon- 
doids are not considered to be good candidates for hy- 
drogen storage applications. However, experimentally it 
has been shown that by treating the surfaces of diamon- 
doids with chemical solutions 3 ^—, the hydrogen atoms 
of diamondoids can be selectively replaced by other com- 
pounds or elements, resulting in significant changes in 
their electronic structure o 37 ! 38 . Among the experimen- 
tally formed diamondoid-based complexes, one modified 
with Li/Li +17 i 39 ' 40 can attract attention for hydrogen 
storage applications and it would be worthwhile to con- 
sider their storage properties in detail. It should be noted 
that it has already been shown theoretically that Li can 
improve the hydrogen storage property of carbon materi- 
al a 41 ' 45 ~ — . However, experimentally it is still controver- 
sial. The wide range of experimental results for hydrogen 
storage improvements after Li doping, from less than 1 
wt% to several tens of weight percentages at moderate 
pressures and temperatures d2rdi ; are mainly attributed 
to sample preparation issues 4 ^—. 

Since ADM has four acidic hydrogen atoms, it can 
possibly be modified with four Li atoms/ions. Our cal- 
culations show that in ADM molecules, more than one 
acidic hydrogen cannot be replaced by a Li + because the 
structure of ADM is significantly deformed such that the 
binding energy of the second substituted Li + becomes 
positive, +0.15 eV . Figure Q] shows the most stable 
configurations of ADM when one of its acidic hydrogen 
atoms is replaced by Li/Li + . From Fig. [TJ it is observed 
that in ADM. Li, Li forms a bond with the host carbon 
atom, similar to other acidic hydrogen atoms, while in 
ADM.Li+, the Li+ is tilted toward the neighboring hy- 
drogen atoms. There is another possible configuration 
for ADM. Li (ADM.Li+) with a tilted (straight) Li-C 
(Li + — C) bond, which is ^0.05 eV less stable than the 
one shown in Fig. [1] Therefore, it is expected that both 
configurations with straight and tilted Li/Li + — C bonds 
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FIG. 1: (Color online) DFT-optimized structures of (a) 
ADM. Li, (b) ADM.Li + and their excess (red) -depletion 
(blue) charges and the related total and projected densities 
of states. 



can be formed in ambient experimental conditions. Ta- 
ble U shows the consecutive binding energy of Li/Li + in 
ADM.Li m (m=l-4) and ADM.Li+. The trend of the 
changes of binding energies of Li to ADM.Li m are ob- 
served to be the same using different methods, PW91, 
PBE, M05-2X, and MP2; by increasing the number of Li 
atoms, the binding energy of Li to ADM.Li m decreases. 
Furthermore, our calculations show that the binding en- 
ergy of Li in ADM. Li is almost same as the binding en- 
ergy of Li + in the ADM.Li + complex. To confirm the 
stability of ADM.Li m /Li + , we have calculated the vi- 
brational spectra of our designed structures. The ob- 
tained results show the absence of any imaginary fre- 
quency, indicating that the optimized structures are real 
minima. We also performed MD simulations at 400 K 
for 48 ps and observed no detachment of Li/Li + from 
the structures and no Li aggregation on an individual 
ADM.Li m complex. As previously mentioned, Li/Li + has 
two nearly isomer configurations with straight and tilted 
Li/Li + — C bonds; as a result, during the MD simulations 
at high temperature, Li/Li + has pendulum movement 
between two local minima. The possibility of formation 
of ADM.Na and ADM.K complexes was also considered, 
and the binding energy of Li to ADM (see Table U) was 
found to be larger than that of Na (-0.86 eV) or K (-0.76 
eV) atoms to ADM. Hence, we focused our examination 
on just ADM.Li m /Li + complexes. 
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TABLE I: E b , the calculated consecutive binding energies of 
Li /Li + in ADM.Li m /Li + complexes. Single point MP2 calcu- 
lations were done for the optimized structures obtained from 
M05-2X/6-311+G(2df,p). 



Cluster 


PW91° 


PBE" 


M05-2X 6 


M05-2X C 


MP2 d 


E b 




E b a 


E b 


dhi-C 


E b 


E b 


ADM. Li 


-1.46 


2.028 


-1.29 


-1.57 


2.010 


-1.51 


-1.59 


ADM.Li 2 


-1.16 


-2.050 


-1.02 


-1.36 


2.023 


-1.30 


-1.36 


ADM.Li 3 


-0.94 


2.075 


-0.84 


-1.17 


2.041 


-1.14 


-1.24 


ADM.L14 


-0.71 


2.094 


-0.62 


-1.06 


2.059 


-1.05 


-1.25 


ADM.Li + 


-1.62 


2.082 


-1.60 


-1.52 


2.038 


-1.46 


-1.43 



"Basis set: plan-wave. 

'DFT cluster method, basis set: 6-311+G(2df, p). 

C DFT cluster method, basis set: 6-31+G(d, p). 
6 Basis set: 6-311+G(2df, p). 



To consider the nature of Li/Li + bonding in 
ADM.Li/Li + , we performed the excess and depletion 
charge analysis and plotted their total and projected den- 
sity of states, as shown in Fig. [TJ The projected density 
of states was plotted for the Li/Li + and the host car- 
bon atoms. From excess and depletion charge analysis, 
in ADM. Li, Li was found to become positively charged 
by donating its 2s electron to the cluster, particularly the 
host carbon atom and all the hydrogen atoms in the com- 
plex. In ADM.Li + , Li + becomes less positively charged 
by accepting electrons, mainly from the hydrogen atoms 
of the cluster. From the projected density of state cal- 
culations, the electrons donated to Li + were observed to 
move not only to its 2s orbital, but also to its 2p y orbital. 
In the ADM. Li complex, the 2s orbital of Li is partially 
hybridized with the 2p z orbital of the host carbon atom 
(see PDOS peaks at -0.25 eV) and makes a sp 3 -like bond 
with it. In ADM.Li + , the 2s and 2p y orbitals of Li + are 
partially hybridized with the 2p z orbital of the host car- 
bon atom (see the states close to Fermi energy). From 
our calculations, it is concluded that the bonding nature 
of Li/Li + in ADM.Li/Li + is predominantly ionic and par- 
tially covalent. 

Before considering the hydrogen storage properties of 
ADM.Li/Li+ complexes, let us consider the possibility of 
clustering of these new structures when two or more of 
them get close to each other. The ADM. Li complexes can 
form dimers without any energy barrier when they get 
close to each other along their heads with Li. From our 
projected density of states analysis, when two ADM. Li 
approach head-to- head (see Fig. [Ha)), it was observed 
that the 2p y orbitals of two Li atoms are hybridized 
with a binding energy of -1.95 eV. Our calculations show 
that ADM. Li complexes can also become connected as 
a chain-like structure with a binding energy of -0.33 eV, 
see Fig. [2^b). When the Li at the head of one of the 
ADM. Li gets close to the ending hydrogen atoms of the 
other complex, they are bonded due to their electrostatic 
interaction between the positively charged Li of the first 
complex and the negatively charged hydrogen atoms of 
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(b) 



FIG. 2: (Color online) DFT-optimized structures of (a) dimer 
and (b) chain-like configurations of two attached ADM. Li 
complexes. 



the second complex; see the excess-depletion charges in 
Fig. [TJ For the chain formation, there is also no energy 
barrier. However, our energy calculations show that the 
dimer configuration is much more stable than the chain 
structure. Therefore, it is expected that dimers would be 
more abundant than the chains in the experiment. In the 
case of ADM.Li + , no clustering was observed. This is be- 
cause in this type of complex, either Li or hydrogen atoms 
are positively charged. Hence, there is a strong repul- 
sion between two ADM.Li + complexes that keeps them 
away from each other and prevents their clustering. It is 
worth mentioning that Xue and Mansoori have recently 
found that ADM.Na + complexes are self-assembled like 
a molecular crystal^ by performing MD simulations on 
125 ADM.Na 4 " complexes. The vacant spaces between 
the complexes may make it possible to store hydrogen 
molecules in high gravimetric weight percentages. More- 
over, recently, experimentalists have been able to prepare 
the positively charged alkali metal doped MOF systems 
via electrochemical reductio n 54 i 55 . They proved that the 
hydrogen storage capacity increases after introducing the 
alkali metal charge cations. Our studied system is an or- 
ganic molecule, which can be easily charged electrochem- 
ically than MOFs. Hence we believe that ADM.Li + is a 
superior to ADM. Li as a candidate for hydrogen storage. 
However, we would like to continue evaluation of both 
ADM. Li and Li + structures, because recent advances in 
modifications of ADM in experiments may make it possi- 
ble to separate ADM.LL„ nano-particles from each other. 

The next step of our study was to consider the hy- 
drogen storage property of ADM.Li + and ADM.Li m 
complexes. As summarized in Table [TT] and as seen 




FIG. 3: (Color online) DFT-optimized structures of (a) 
ADM. Li and (b) ADM.Li+ when one to five hydrogen 
molecules are adsorbed. 



from Fig. [3J each Li/Li + site adsorbs a maximum of 
five H2 molecules. Therefore, it can be predicted that 
the gravimetric weight percentage of hydrogen storage 
for ADM.Li+ is -7.0 % and between 7.0-20.0 % for 
ADM.Li m complexes, if experimentalists can find a way 
to prevent their clustering. If clustering occurs the stor- 
age properties of ADM.Li m will be less than the above 
mentioned value. 

Now let us focus on the hydrogen adsorption on 
Li/Li + functionalized ADM. The calculated binding en- 
ergies along with the bond parameters for the hydrogen 
molecules using various functionals and basis sets are 
given in Table [TTJ From the table, it is seen that in 
our case study, both pure and hybrid functionals pro- 
vide binding energies close to the more accurate MP2 
method. It is observed that the binding energies of hy- 
drogen molecules on ADM.Li m /Li + are on the order of 
-0.1 to -0.23 eV, which is very good for hydrogen storage 
applications. Furthermore, calculated binding energies 
for the cationic Li + are higher than those of the neutral 
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TABLE II: Eb (in eV), Lh-h, dn 2 -Li, and dLi-c (in A) are the calculated binding energies of adsorbed hydrogen molecules, 
bond length average of hydrogen molecule(s), average distance between the center of hydrogen molecule(s) and Li/Li + , and 
average bond distance of Li and host carbon atoms in different methods, respectively. 
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PBE a 




M05-2X 6 




M05-2X C 


MP2 d 


E b 


- 

Lh-h 


- 

dn 2 -Li 


dhi~C 


E b 


E b 


Lh-h 


- 

dH 2 -Li 


dLi-C 


E b 


E b 


ADM.Li-(H 2 )i 




-0.10 


0.753 


2.184 


2.020 


-0.10 


-0.11 


0.743 


2.149 


2.008 


-0.11 


-0.11 


ADM.Li-(H 2 ) 2 




-0.20 


0.788 


1.788 


2.007 


-0.20 


-0.15 


0.760 


1.868 


1.948 


-0.12 


-0.10 


ADM.Li-(H 2 ) 3 




-0.15 


0. / IU 


1.882 


2.049 


-0.14 


-0.15 


0. (49 


1.900 


2.011 


-0.14 


-0.12 


ADM.Li-(H 2 ) 4 




-0.14 


U. / Do 


l.yoo 


2.088 


-0.13 


-0.15 


n 7/1 q 
U. / 4o 


O flQ/1 
z.Uo4 


2.035 


-0.14 


-0.12 


ADM.Li-(H 2 ) 5 




-0.11 


U. /OU 




2.085 


-0.11 


-0.13 


U. ( 40 


O 1 07 
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-0.13 


-0.10 


ADM.Li 2 -(H 2 ) 
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-0.14 
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-0.13 


-0.14 
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-0.14 


-0.11 


ADM.Li 3 -(H 2 ) 
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-0.14 
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-0.13 


-0.14 
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ADM.Li 4 -(H 2 ) 
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-0.14 
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2.065 
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-0.10 


ADM.Li 4 -(H 2 ) 


20 
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-0.11 


-0.12 
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0.746 
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2.046 


-0.17 
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ADM.Li+-(H 2 ; 
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-0.16 


-0.20 


ADM.Li+-(H 2 ; 
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-0.19 
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2.049 


2.184 


-0.17 


-0.19 


0.746 
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2.145 


-0.15 


-0.19 


ADM.Li+-(H 2 ; 
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-0.17 


0.755 


2.139 


2.223 


-0.15 


-0.17 


0.744 


2.124 


2.187 


-0.14 


-0.17 


ADM.Li+-(H 2 : 


Is 


-0.15 


0.755 


2.236 


2.270 


-0.13 


-0.16 


0.744 


2.192 


2.241 


-0.13 


-0.15 


H 2 






0.749 










0.739 











"Basis set: plan-wave 

'DFT cluster method, basis set: 6-311+G(2df, p). 

C DFT cluster method, basis set: 6-31+G(d, p). 
d Basis set: 6-311+G(2df, p). 



Li-doped ADM system. The calculated Li-H distance 
increases with an increase in the number of hydrogen 
molecules. The small changes in binding energy of hy- 
drogen molecules can be attributed to various reasons: 
the amount of positive charges on Li£i, the distances be- 
tween H2 molecules and Li, the Li— C bond distance^, 
the interaction between the hydrogen molecule a 52 ' 53 , the 
strength of hybridization of hydrogen molecules with Li, 
etc. It is observed that the adsorbed H 2 molecules are 
located at distances of ^2.1 A of Li/Li + . It is also 
seen that the first H2 is adsorbed on top of the Li/Li + 
(see Fig. 13]), but when the second, third, or fourth H2 
molecules are adsorbed, they prefer to move to the lateral 
side of Li/Li + . As observed from Fig. [31 when the hy- 
drogen molecules are adsorbed on the ADM.Li/Li + com- 
plexes, the position of Li/Li + changes between their two 
local minima with a straight or tilted Li/Li + — C bond 
in order to reduce the steric repulsion between the ad- 
sorbed hydrogen molecules and the hydrogen atoms of 
the ADM.Li/Li+ structures. 

To find the nature of bonding between H2 and Li/Li + , 
as an example, we considered the excess and depletion 
charge and the projected density of states for ADM. Li 
when one or four H2 molecules are adsorbed on a Li 
atom (see Fig. [3}. From the excess and depletion charge 
analysis, positively charged Li was shown to polarize H2 
molecules under its induced electric field. Due to this 
polarization, there is a small bond elongation for the H2 
molecules, as seen from Table HH but no dissociation of 



H 2 molecules. Density of states analysis indicated that 
when the number of H2 molecules increases from one to 
four, they start to interact with each other such that the 
states related to H 2 molecules are broadened^, between 
—8.0 eV and -10.0 eV. Regarding the ADM.Li-(H 2 )i, as 
shown in Fig. 01 the Is orbital of H2 molecule is slightly 
hybridized with the p z orbital of the Li atom at -10 eV. 
When the number of H2 molecules changes to four, the Is 
orbitals of H2 molecules prefer to interact with the p x and 
p y orbitals of Li (at energies between ~-8.0 eV and -8.5 
eV) instead of its p z orbital. Hence, when the number of 
H 2 molecules increases, they prefer to locate in the lateral 
positions of Li rather than in the top position. To better 
understand how the electric filed induced by positively 
charged Li affects the binding energies of H 2 molecules, 
we plotted the magnitude of the induced electric field 
along the Li/Li + — C bond for ADM.Li/Li + , shown in 
Fig As shown by the figure, the amount of generated 
electric field at the center of adsorbed H2 molecule on 
Li/Li+ is on the order of 2.1/3.4 (V/A). The polarizabil- 
ity of a hydrogen molecule along (ck||) and perpendicular 
to the hydrogen molecule axis in an external elec- 
tric field are 6.3 a.u. and 4.85 a.u., respectively^— . 
The adsorption energy of hydrogen in such electric fields 
is estimated to be E 6 ~-l/2 P.E ext ~- l/2a £f xt ~-0.11/- 
0.29 eV. These values are very close to ones reported in 
Table HT1 Therefore, it is expected that the electrostatic 
interactions between H2 molecules and Li make larger 
contributions to the binding energy of H2 molecules on 
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(a) ADM.Li-(H 2 ) 1 



(b) ADM.Li-(H2)4 
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FIG. 4: (Color online) Excess (red)- depletion (blue) charge 
iso-surfaces and total and projected densities of states for (a) 
ADM.Li-(H 2 )i and (b) ADM.Li-(H 2 ) 4 . 
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FIG. 5: (Color online) Induced electric field in the direction 
of the Li/Li+-C bond of ADM.Li/Li+. Li/Li+ located at the 
origin. The dotted lines indicate the center of the adsorbed 
hydrogen molecule. 
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ADM.Li/Li + complexes than their hybridizations. It is 
worth mentioning, as seen from Table HH that the aver- 
age of adsorption energy of H 2 molecules on ADM.Li m 
does not change significantly when the number of Li sites 
increases. This indicates that in such complexes, the 
interaction between hydrogen molecules and positively 
charged Li is still electrostatic. It should be mentioned 
that Yoon et al~ have used a similar analysis and suc- 
cessfully explained why hydrogen molecules can be ad- 
sorbed on charged and doped carbon fullerenes. Also, 
Zhou et al£Q have recently shown that under an applied 
electric field, hydrogen storage property of a boron ni- 
tride sheet is substantially improved due to polarization 
of the hydrogen molecules as well as the substrate. 

To investigate the thermodynamics of adsorption of 
H2 molecules on ADM. Li and ADM.Li + , the occupation 
number of H2 molecules per site (Li/Li + atom) was cal- 
culated as a function of the pressure and temperature 
using the following formul a 61 ' 62 : 

= n=o , 2) 

n=0 



FIG. 6: (Color online) Occupation number as a function of the 
pressure and temperature on (a) ADM. Li and (b) ADM.Li + . 



where (i is the chemical potential of the H2 gas, e„(< 
0) and g n are the (average) binding energy of the H2 
molecules and the degeneracy of the configuration for 
a given adsorption number of the H2 molecules n, re- 
spectively and, k and T are the Boltzmann constant 
and the temperature, respectively. Figure |5] shows the 
occupation number of H 2 molecules on the Li and Li + 
atoms as a function of the pressure and temperature 
where the experimental chemical potential of H2 gas 63 
and the calculated binding energy (e„) obtained from 
MP2 calculations were used. The occupation number 
of H2 molecules / at 150 K and 30 atm in both cases is 
~4, as shown in Figure [6l This is attributed to the Gibbs 
factor (e 4 ( /i ~ C4 )/ fcT ) for the binding of four H2 molecule, 
which dominates at 150 K and 30 atm (/j = —0.08 
eV, £4 on the Li and Li + is -0.12 and -0.17 eV, respec- 
tively). The number goes to zero at room temperature 
(/j =n*> —0.32 eV). Therefore, this analysis shows that the 
ADM.Li/Li + structures may have considerable potential 
as high-capacity hydrogen storage materials. 
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IV. CONCLUSIONS 

Using different first-principles approaches, we have 
shown the possibility of formation of ADM. Li and 
ADM.Li+. It is predicted that ADM.Li+ structures will 
not be clustered under ambient conditions, while there 
is a tendency for clustering of ADM. Li complexes. We 
found that each Li + is capable of holding five hydrogen 
molecules in molecular form. Binding energy of the hy- 
drogen molecules was in the range of -0.15 eV to -0.23 
eV, which was due to the induced electric field created 
by the Li + ion. Furthermore, the calculated occupancy 
number reveals that ADM.Li/Li + to be ideal candidates 
for hydrogen storage. Thus, ADM.Li + is expected to be a 



superior candidate for use as a hydrogen storage medium, 
as hydrogen molecules can be adsorbed or desorbed from 
this complex at pressures and temperatures accessible in 
recent technologies. 
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